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CHAPTER  I 


INTRODUCTION 
I.  Introduction 

This  dissertation  presents  three  examples  of  the  use  of  a 
collision  term  derived  through  the  Klimontovich  formalism  for  the 
calculation  of  transport  coefficients  of  a plasma  driven  by  a high- 
frequency  electric  field  which  is  not  necessarily  weak.  The  resistive 
heating,  which  has  been  calculated  before,^  is  derived  here  using  the 
formalism  of  Goldman  (1970),  as  is  the  electron-ion  thermal  energy 
exchange  rate.  The  collision  term  as  given  by  Klimontovich  and  Puchkin 
(1974)  is  used  as  the  basis  of  an  expression  for  the  thermal  conduc- 
tivity. The  presentation  is  simplified,  in  most  places,  by  supposing 
the  ions  to  be  singly  charged.  In  these  examples,  the  effect  of  an 
ambient  magnetic  field  is  not  considered. 

Chapter  I introduces  some  topics  which  will  be  used  in  the 
computations.  Chapter  II  outlines  the  derivation  of  the  collision 
term  used  in  the  computations . The  computation  of  an  inverse 
Bremsstrahlung  rate  and  a heat-exchange  rate  from  the  frequency-domain 
representation  of  Goldman  is  done  in  Chapter  III.  The  electron 
thermal  conductivity  is  computed  with  the  use  of  the  Chapman-Enskog 

*More  complete  references  for  these  topics  are  given  in  the 
specialized  chapters. 
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method  in  Chapter  IV.  The  results  of  the  computations  of  Chapters  III 
and  IV  are  Eqs.  (III. 41),  (III. 64),  and  (IV. 71).  Their  values  are 
shown  in  Figures  3,  4,  and  5 and  in  Tables  I and  II;  the  most 
important  result  is  probably  the  expression  for  the  thermal 
conductivity. 

Gaussian  c.g.s.  units  are  used.  Use  of  the  symbol  y'.  does  not 
imply  y an  integer,  and  —I,.  is  defined  by  continuity  if  yl  doesn't 

y • 

exist.  Many  symbols  in  the  text  have  meanings  which  are  conventional 
or  inferrable  from  the  context,  and  are  not  defined  when  introduced. 

| Most  of  the  symbols  used  are  defined  in  Chapter  VI. 

2.  Transport  Processes  in  a Plasma 
A thermodynamically  isolated  plasma  not  in  equilibrium  will  tend 
to  move  toward  equilibrium.  If  the  departure  of  the  system  from  an 
equilibrium  configuration  is  small  enough,  and  the  equilibrium 
configuration  is  not  a special  point  for  the  perturbing  or  restoring 
agents,  then  the  rate  of  restoration  will  be  proportional  to  the 
perturbation.  If  an  infinitesimal  perturbation  is  maintained  by 
interaction  of  the  plasma  with  outside  forces,  the  plasma  will  reach 
a stationary  state  in  which  the  perturbation  is  balanced  against  the 
plasma  response.  The  constant  of  proportionality  between  the 
departure  from  equilibrium  and  either  the  perturbing  force  or  the 
rate  of  restoration  is  called  a "transport  coefficient." 

The  behavior  of  systems  which  deviate  infinitesimally  from 
equilibrium  is  discussed  in  the  books  by  Prigogine  (1967), 
deGroot  (1951),  and  Cox  (1955). 
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Examples  of  departure  from  equilibrium  are  nonthermal  distri- 
bution functions,  temperature  differences  between  species,  macroscopic 
electric  fields,  gradients  of  temperature  or  concentration.  The 
restoration  rate  is  controlled  by  interparticle  interactions,  called 
"collisions"  by  extension  from  the  classical  billiard-ball  molecules. 
Examples  of  transport  coefficients  are  electrical  conductivity, 
interspecies  equilibration  rate,  thermal  conductivity,  rate  of 
relaxation  of  a single  species  toward  thermal,  and  diffusivity. 

As  an  example  of  these  ideas,  suppose  that  two  populations 
of  greatly  differing  mass-say  ions  and  electrons— have  different 
temperatures.  The  coupling  between  species  is  weak  enough  that  each 
population  remains  essentially  thermal,  and  the  temperatures  will 
tend  toward  equality  according  to 


The  olectron-ion  equilibration  time  is  • transport  coefficient 
for  the  process.  This  Is  an  example  of  a "transport"  process  which 
isn't  a flux  in  configuration  space. 

Another  example  of  a system  with  transport  is  a plasma  with 
imposed  temperature  gradient  and  electric  field.  A stationary  state 
is  described  hy 

J = €rE  + o<  V9  (2a) 


£ =-/£  - KV&  (2b) 

giving  the  electrical  current  J and  the  heat  flux  S (Spitser  and 
Harm  1953).  The  transport  processes  here  are  the  actual  fluxes  of 
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{ heat  and  charge,  and  the  transport  coefficients  are  €" , «* , O,  and 

I K 

I 

In  a system  with  a finite  departure  from  equilibrium,  the  ratio 
of  rate  of  restoration  to  strength  of  perturbation  is  likely  to  be 
positive,  though  perhaps  not  constant.  If  a periodic  electric  field 
is  applied  to  a plasma  with  species  of  different  temperatures,  the 
species  will  still  exchange  thermal  energy,  though  the  process  is 
no  longer  one  of  equilibration.  If  the  field  is  strong,  it  can  modify 
the  response  to  an  applied  (constant)  electric  field  or  temperature 
gradient;  the  thermal  and  electrical  conductivity  are  still  constant 
with  respect  to  the  perturbations,  but  are  functions  of  the  driving 
field.  If  the  plasma  response  considered  is  that  of  the  plasma  to 
the  driving  field  itself--the  rate  of  heating  by  inverse  Brerasstrahlung- 
then  the  response  is  nonlinear  when  the  driving  field  becomes  strong. 

3.  The  Effect  of  a Driving  Field 

When  a plasma  is  so  strongly  driven  by  an  electromagnetic  field 
that  the  resultant  electron  velocity  much  exceeds  the  thermal  speed, 
the  most  important  effect  may  be  the  generation  of  plasma  instabilities 
(DuBois  and  Goldman  1965;  Silin  1965)  and  electrostatic  turbulence 
(Kadomtsev  1964;  Tsytovich  1972).  But  the  turbulence  needs  time  to 
grow,  there  may  be  effects  which  are  independent  of  turbulence,  and 
it  is  possible  for  the  driven  speed  to  be  near  the  thermal  speed  with- 
out turbulence  (Brownell,  Dreicer,  Ellis,  and  Ingraham  1974  and  1975). 

An  alternating  electric  field  will  change  the  average  relative 
speed  of  the  two  species  and  thereby  influence  transport  processes. 

The  effect  will  become  Important  as  the  driven  speed  of  the  electron 
approaches  its  thermal  speed.  A good  discussion  of  strong-field  effects 
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is  given  by  Kidder  (1974). 

The  intensity  of  the  electromagnetic  wave  which  drives  the 
electrons  at  their  thermal  speed  depends  on  the  electron  temperature 
and  the  wave  frequency.  It  is  shown  in  Figure  1.  For  that  figure, 
it  is  supposed  that  the  wave  is  propagating  at  the  vacuum  speed  of 
light;  physically,  it  could  be  important  that  the  field  strength 
increases  in  a wave  approaching  a reflection  point.  The  intensity  at 
which  the  nominal  driven  speed  reaches  the  speed  of  light,  so  that  the 
treatment  given  here  needs  modification,  is  shown  in  Figure  2. 

4.  Integration  with  respect  to  Wave  Number 
The  expressions  derived  in  Chapters  III  and  IV  involve  velocity 
moments  of  the  collision  term  (11.29)  in  the  next  chapter,  and  those 
moments  involve  integrals  over  k-space  (and  sums  over  the  integers  n) 
of  expressions  weighted  by  exp  For  large  values  of  k, 

the  integrands  are  proportional  to  k . We  can  expand  each  integrand 
in  a series 


I = 4*f  I.*  I, 


(3) 


' and  take  the  leading  term.  In  this  process,  the  cutoff  has  disappeared 

i 

that  was  imposed  by  the  exponential  at  small  values  of  k,  so  we  re- 

I 

i impose  it : 


Jw1  - 

The  expression  on  the  end  contains  the  Coulomb  logarithm 
which  involves  the  limits  on  k.  The  lower  limit  k,,,^  will  be  an 


(4) 


(5) 
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typically  at  least  5 or  10  (Spitzer  1962,  p.  128),  the  logarithm  can 
be  estimated  by  ignoring  small  multipliers  in  k^j^,  and  we  can  use 


L ~ L 
M*»/n  *0  Up  • 


(6) 


The  upper  limit  k,,,^  is  the  inverse  of  the  distance  of  closest 
approach  b0,  and  is  required  by  some  terms  ignored  in  the  derivation 
of  the  collision  term  in  Chapter  II.  (A  brief  discussion  is  given  in 
Appendix  F.) 

Another  form  for  (4)  is 

3 


/<V" 


fy>v~e 


He*Z  2(2  if)"1 


jV-- 


(7) 


ua  is  the  conventional  electron-ion  collision  frequency  for  momentum 


transfer  (Perel  and  Eliashberg  1961;  DuBois,  Gilinsky,  and  Kivelson 
1963;  Kidder  1971): 

nyfzTf  n;  e*e*  /CrvA 
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(8) 


What  happens  to  A in  a strong  alternating  field?  The  minimum 
impact  parameter  decreases  as  v"2,  and  the  effective  quantum-mechanical 
size  of  the  electron  only  as  v"1,  so  that  at  some  energy  the  effective 
distance  of  closest  approach  is  determined  by  the  electron  size  instead 
of  b0.  A parallel  argument  is  siven  by  Marshak  (1941).  Moreover,  the 
effective  speed  is  some  average  of  the  thermal  speed  and  the  driven 
speed;  so  knax  can  be  calculated  by 


. 


I 
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'majr 


m"a-‘  K .81’ 


V max  (v„  ) 

Vpp  1 ff*  mu/  ' 


If 
! • 


t = 


fmv 


(9) 


JeB  toVeFF  ~ i^^FF 

The  crossover  to  quantum-mechanical  limitation  happens  when  2«.c/ve££, 
in  which  OC  is  the  fine-structure  constant,  decreases  past  1;  and 
that  happens  when  the  electron  energy  increases  past  about  40  ev. 

The  effect  of  the  strong  field  on  k,^  can  be  ignored  within 
our  crude  approximation,  because  although  the  effective  shielding 
distance  may  be  increased  by  decorrelation  of  the  ‘.wo  species,  each 
species  is  still  shielding  itself. 

In  a plasma  with  turbulence,  Chere  is  often  important  behavior 
in  the  spectrum  of  (e.g.,  electron  density)  fluctuations  corresponding 
to  propagating  Langmuir  waves  (see,  e.g.,  DuBois  and  Goldman  1972), 

and  that  spectrum  is  strongly  damped  unless  k£0.3kD.  Thus  weak 

/*o 

dk,  and 

r 

non-coleective  collisions  contribute  / dk,  where  an  upper  limit  may 
have  to  be  imposed.  A variation  of  this  idea  is  discussed  by  Kihara 
and  Aono  (1963). 

5.  The  Quantities  f^T(^) 

This  section  defines  and  gives  some  properties  of  a family  of 
functions  that  enter  the  analysis,  and  which  have  been  given  the 
labe 1 Pr . 

Definition; 

TO)  = Jim  E Mf  *■  * 

■7-*°  n s -co  * / 
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Expression  for  T r: 

A derivation  of  an  expression  for  P is  given  here.  (Another 
derivation  is  given  in  Appendix  A.)  The  starting  point  is  the 
definition  of  T0: 

pa)  = z:  an 

0 7 7T  " *7 

An  easily  derived  identity  (Watson  1944,  p.  32)  is 
. nV/Z 

= J.£  J'dSfr*-'- >”’*<*  (ID 

For  (fixed)  nonzero  the  contribution  from  the  JQ  term  will 
come  only  from  arguments 

B'/'^  ± &(l).  (13) 


When  *1  is  small  enough,  the  contributing  part  is  a smooth  function 
of  n,  and  the  sum  is  tractable: 

* s nin>  £ ~^I/r  _ <u) 

rji)  = (15> 

in  which  an  integral  representation  for  JQ  is  used. 

J Jzc~xtcos  = TT^t0^  (16) 

c«  - iS&r'Jk  **'-*/*<»#•**  (id 

0 / • 

PQ  is  real,  so  2 JdO  -> Jd 0. 

• “•  « 

PW  = f ^ 7T*  C*f>  [-(&*}*  2i  Woe]  (I8) 

0 •'a  7**  1* 


in  which  a modified  Bessel  function  appears: 

= e'*Ih(x) 


(21) 


An  expression  for  the  other  P's  can  be  obtained  by  including  a 
factor  p in  the  argument  of  the  exponential  in  (11).  The  value  of 

. 2i 


the  expression  becomes  ) » 80  thilt 

= f-rpApYi8; 

From  Abramovitz  and  Stegun  (1964,  p.  377),  we  have 

A (**)  = 


(22) 

(23) 


(24) 


and  also  that 
fh) 


(n't)! 


\am  « c-t)  m -«•;  <25> 

to  give  the  result 

The  symbol  M is  the  confluent  hypergeometric  function  (Jahnke  and  Emde 
1945,  p.  275). 


i 


11 


A related  quantity: 

The  expression  preceding  the  equality  *ign  in  (33)  below  enters 
the  analysis.  It  can  be  signified  by  the  change  of  variables 

7 ■ <27: 

i = yc£  <28: 

in  which  a nondimensionalized  field  strength  is  introduced: 

E = — = - . (29; 

ZYC 

The  argument  of  the  Bessel  functions  is 

k’slt  = kyu  ~ = . (3o; 

The  susceptibility  term  is 

I1  /"tf* 

A*  e‘  o« 


, iSt 


IT  * 


Now  the  expression  of  Interest: 

ZL  ^ Av-./i"  ) Qvn,yc(n<*>o)  Jy,  (k’dt') 

O c 

= Zfr’^^c-T’Xaji.)  <33> 

Operating  on  that  expression  with  k-2  jlim^k2  gives  the  result 

® 2.^1 

k-tID  n*-<D  * — 

« r-r)V'*£m  <34) 
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CHAPTER  II 


THE  COLLISION  TERM 


The  evolution  in  time  of  a particle  distribution  function  obeys 


the  Boltzmann  equation 


iL  +y.  i£ 


The  right-hand  side,  the  "collision  term,"  is  that  part  of  — 
which  is  generated  by  interparticle  interactions.  In  this  chapter, 
expressions  are  developed  following  the  treatment  of  Goldman  (1970) 
for  two  velocity  moments  of  the  collision  term 

<: 

in  which  the  Fourier-transform  variables  k and  cj  appear.  An 
expression  given  by  Klimontovich  and  Punchkin  (1974)  for  the  time- 
dependent  form  is  noted  for  later  use. 

1.  Development  in  the  Time  Domain 
A specified  set  of  identical  point  particles  can  be  described 
by  the  exact  ("Klimontovich")  distribution  function 

Mie  " 

f <%&)  - L « 

/■»/ 

where  (^(tj.nrv^t))  is  the  position  of  tne  Ith  of  the  N particles 
in  a (six-dimensional)  single-particle  phase  space.  It  has  the 


l 


14 


exact  equation  of  motion 

( +x-jg)P  (z,v;t)  +£-•[»» (firwof  J = 

This  equation  is  useful  if  the  right-hand  side  is  set  to  zero. 

&mc  (x,v,t)  is  the  acceleration  that  the  field  at  (x,v,t)  would 
give  a particle  there;  the  self-field  of  the  particle  is  to  be 
neglected.  A parallel  development  using  the  wave  equation  and  a 
semiclassical  treatment  of  the  field  seems  feasible  (Rand  1964;  Brown 
and  Kibble  1964;  Pert  1972;  Seely  and  Harris  1973;  Seely  1974),  but 
is  more  cumbersome  in  describing  classical  particles. 

The  system  can  be  considered  to  be  a member  of  a Gibbs  micro- 
canonical  ensemble,  and  averages  taken  over  the  ensemble.  We  will 
use  the  notation  suggested  by 

P(x,yri)  = < p”"(x,Xt)>,  (5a) 

= P “f.  (5b> 

Taking  the  ensemble  average  of  each  term  of  (4)  with  ti=0  leads  to 

(fr  = -h"<s&r*rl)sP(irt)>>  (6) 

and  a problem  in  kinetic  theory  is  the  evaluation  of  the  right-hand 
side.  The  Boltzmann  collision  integral  and  that  of  Landau  are  based 
on  the  picture  of  a binary  collision,  and  a statistical  argument  leads 
to  the  Fokker-Planck  form  for  (6)  (Chapman  and  Cowling  1939,  p.  46  ff; 
Chandrasekhar  1943;  Wu  1966,  p.  32;  Present  1958;  Landau  1936;  Krall 
and  Trivelpiece  1973;  Delcroix  1960,  p.  116). 


L5 


The  difference  of  (4)  and  (6)  is 

(i-  +r-2-)tfrzri>  +j?-!isP  = -fr-[rsa.+sfsg-<tdt>]y) 

9 1 " 

We  specialize  now  to  a gas  of  particles  coupled  only  by  electromagnetic 

g 

forces;  then  commutes  with  . We  ignore  terms  quadratic  in  the 
fluctuations.  Allowing  more  than  one  kind  of  particle  in  the  system 
and  introducing  a species  subscript  * gives  the  equations 

(k  * JV’  <8> 

c & +£<&)£(&)  * y««  h <*> 

The  treatment  in  this  chapter  follows  that  of  Goldman  (1970). 

A Green's  function  for  the  operator  in  (8)  is 

&(t-t')Qjzzvr'it')  - 8a-m[z (10) 

in  which  x°  and  v°  are  the  coordinates  a particle  would  have  at 
time  t if  it  began  at  (x/v't*)  and  moved  with  acceleration  0*  (x,v, t): 

Y*(t;zr V)  = y'  + f <4r  a (z (ri  vfr)  r)  dD 

n*  r T> 

r'(t;zYt')-  x'+jr'a-iO  WW  £ (z(rt)tfTt)  (12) 

•f  ' y * 

G„has  the  property 

<oj:*iYV  tr)  - Sfy-Tr'J'  (13) 

The  formal  solution  of  (8)  which  takes  a specified  value  at  t « tQ 


s Part)  = \'jVj\'iP(rrt)Q(ii'jt'u) 

of  J 4 o ^ a y 

+ fj’z'jVpt  'm~c)  9 (i'rt')G  H') 

° of  o< 

o 

This  is  not  an  explicit  solution  because  the  source  term  S ^ 
depends  on  through  The  dependence  is  expressed  by 


Poisson's  equation 


V-£F(x<)  = = WH  s?(xi) 

«=$/  < 

= «rlf,  py  sfxzyt).  (1 

If  we  define  a homogeneous  plasma  as  one  for  which  the  ensemble 
averages  ? (rpt  Jand  ff/y»/)are  independent  of  and  suppose  our 
system  to  be  such  a plasma,  then  C depends  on  X and  Z only 
through  the  difference  Z‘%',  and  we  can  take  space  Fourier 
transforms: 

sf/vij  5 pGe-‘^x  sP(iv  i)  ci 

6 (ft'il'J  = a 

Offc  ' J ~ « 

The  quantity  SE^O)  has  an  important  property:  ££(*,0  is  the 
field  of  a set  of  point  particles  moving  at  nonrelativistic 
speeds,  so  that  interaction  is  effectively  instantaneous  I 
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(19) 


(20) 

(21) 


and  that  la  parallel  to  k,  so  that 

S£k(t)  = %SEa  (t). 

This  defines  the  signed  quantity 
We  use  the  abbreviations 

= t.  [SCt-t-)  + prAuri'] 

V 

Some  Fourier  transforms  are  listed: 

G^Jyv'ti'J  = t(x-x 

sQyij  s fSv'rP/*VGt/t*0 
G^rv'/l'J  - S(v-t'J 

Vy<,“  =W;r'^<27) 

ST  6U  = rf (k  ? o)  (28) 

It  can  be  seen  from  (18)  that  vanishes  at  ^ 

In  terms  of  k-dependent  variables,  the  collision  term  from  (9)  is 


(22) 

(23) 

(21*) 

(25) 


(29) 
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2.  Symbolic  Solution  for  Correlation  Functions 
In  this  section,  expressions  for  ^ and 
for  ($e(0  are  written  down. 

Following  are  some  definitions  for  some  symDols  used. 

The  nonlinear  susceptibility  ^(tjt’): 

-lLy^(U')  = --jk  J/vJv'&t-i'JGiiWtOk'  -^J( 30) 

The  charge  density  fluctuation  for  one  species: 

tyo  = JjV/.  sL^t}  <»> 

The  part  of  the  charge  density  fluctuation  which  has  evolved  from 
the  initial  value: 

H V*  £*(/**>  <32> 

Total  susceptibility  and  charge  density: 

^ 5£i.+*/t  +sp,  <33) 

k <k 9 t rk  'L 

The  dielectric  function: 

ekan  = + \ (tf)  (3i.) 

The  inverse  dielectric  function: 

Jji*€£(tt*)€k(t*t')  - S(t-t')  °5) 

The  last  equation  defines  ; the  questions  of  existence 

k 

and  uniqueness  are  deferred. 

Development  of  the  expressions  for  the  defined  quantities: 


Using  (26)  in  (31)  gives 

V<)  • K(i>  - (36) 

uy.ua  ■-  $-•  6><t-i-)[SvjVG,(K'‘ni  on 

«h  n«  J *<k  )v 

, ?£*  &(t-v) Jj\  ?•  ? ) <38> 

If  tQ  is  chosen  to  be  recent  enough  that  collisions  haven't  had  time 
to  change  the  distribution  function  very  much--i.e.,  if  the  right-hand 
side  of  (6)  is  set  to  zero--then  £ is  explicitly  known  from  its 
initial  value: 

f,m  = r<  t-AKn,,i.) 

If  the  initial  distribution  was  Maxwellian,  this  condition  may  be 
applicable  even  if  collisions  have  had  time  to  be  effective,  because 
the  fastest  effect  of  collisions  is  intraspecies  thermalization. 

Then  using  £(V  - Atij  t‘)  - £(>-^4,0  in  <38>  8ives 

iW  JW)  = tt*  J/y  £•  toi^\ko) 

The  expression  for  E,  from  (28),  is 

t£k<o  - 7rLrr/t)  (41) 

which  is  the  same  as 

Jjteftt’)  SE^t')  = j£-  <42> 


i 
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so  that  an  explicit  expression  for  fE  is 

1^(0-  Jji'e''(tc)  If  (ff(n.  w> 

Me  assume  that  (vtQ)  is  known.  Then  we  can  use  the  expression 
derived  in  Appendix  D following  Wu  (1967): 

*8  * 

This  quantity  is  independent  of  k.  In  Appendix  E,  we  prove  the 
identity 

jl/V'e  (Tcx'ttj G (n'tfl  <45> 

For  the  field  correlation  function,  we  develop 

- J j’v'jV  %)  Gftri l. 1 il- ) 

= /wV  f £ r?-V<,  G,/??'/  4, ) G./W.J 

= 5 ? (r-Ajt.  l0)  GL(i;irtt).  (46) 

«0  * * ' " vk  • - i 

Then 

1 ?.  £ to-  ^ i.  •) f '“r't  ry<' (47) 
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< sr[(t)  f £ jy^r  fcrtje’'**^”***^*'*^ 


This  expression  into  (40)  gives 


<sFt<u  tskcn> 

= Jji -Ji’^an  <;;«  n 

The  term  (fF<P^ca n be  expressed  in  terms  of  this.  Making  the 


substitutions  suggested  by 


jf  = e~‘rr]  ft*  ’ srfp  = + 


SfrfF 


gives 


&P+  <ir/i)  i-f/iv)  (50) 

. % r.  t, ) e-‘f>i(v- 1 

* -t  ) i k-  -i  tr-4i<i  *‘J<fE  (uir  (fl) 

Jx  9 V * *4 

~ (51) 

3.  Specialization  of  the  Driving  Field 
For  the  examples  given  in  Chapters  III  and  IV,  we  specialize  to 
a system  in  a linearly  polarized  electric  field  with  no  magnetic  field: 

*va,of»  - (52> 

It  Is  turned  on  after  time  t in  such  a way  that  for  t>0, 


\it,  » fjrg^tt; t)  = -X 


' cqs  <Jmi 


j 


I 


S'Ji'i/t.  = -4 


In  this,  JJ  « MjL 4 = MlS  . 

•(  r"*  “V  ^ 

The  forms  taken  by  (20)  and  (21)  are 


A^*'  = A*  it,  ~ \ i't*  = -J/0v  ( coso>0t  - (55) 

f/H*'  - 

- -c/^  (:invti-Tmu,ti'J  4 (i-l')U,€  CosW't'  (56) 

= k- [(*;'*■  4u<»*%y)(t't')'e/J[*'"ot  -sin^t'^(57) 


The  susceptibility  (40)  is 


UYjif)  - 


4.  Fourier-Transformed  Expressions 
This  section  gives  a sketch  of  the  derivation  of  two  moments 
of  the  collision  term  for  later  use.  If  the  system  is  stationary 
except  for  the  periodicity  imposed  by  (^(v . t ) periodic  with  period 
eJo,  then  any  function  with  two  time  arguments  f(t,t*)  determined 
by  the  system  can  be  expressed  in  the  form  g(t-t',t'),  where  g is 
periodic  in  the  second  argument.  A double  Fourier  expansion  can  be 


made : 


r<*P  - f 

A* -to  v 


i . 

i i 

li 
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f ! 
• I 

: t 


i1 


fc^-J  -- 

ifdi  - z 


ifUcJ)  = Jji  !f(i)  e 

The  space  and  time  Fourier  transforms  are  normalized  in  the  usual 
way  by  setting  the  (infinite)  intervals  of  integration  to  unity 
when  they  appear  as  factors;  e.g.,  a complete  statement  of  (62) 
could  have  the  form 


(60) 

(61) 

(62) 


r 

e 


- I 
2 


ZTT 


where  the  values  of  T are  integer  multiples  of  — , the  system 

has  period  T,  and  lim  is  implied. 

T 

Because  our  system  may  be  in  fact  not  stationary  but  have 
a time  variation— e.g. , heating  caused  by  the  collision  term— 
which  is  slow  compared  with  the  relaxation  times  of  the  system, 
the  frequency-dependent  coefficients  in  Eqs.  (59)  through  (62) 
may  have  a residual  time  dependence. 

Important  properties  are 

C fa  '{ft*  1 £ & W]  tnM\rw)  = Z ) 

[ pi 'jai') ifdljtn*. ) * Z f (6A) 


1 


25 


The  same  symbol  is  used  for  the  linear  and  nonlinear 

susceptibilities;  they  are  distinguished  by  the  number  of  arguments. 

Making  appropriate  substitutions  leads  to  the  following 
expressions,  which  will  be  used  in  Chapter  III: 


<%.«>  S?k (t)>  Uv.) 

-'tiLjgJ  f<{ *, 

/r?/-  * L to  i 


<i  <j/>  ^.5/ 

xti-'(*k)*jn-J,)uk-i)  (72) 

ny»»  *-/•/>}  fjj  J 


* - “I 

X z Mr„  ( i “)  (k^f 

+ ^r-  M-  [(*  •i*~h  ) 

* £ Mr„ (k- % .<W  ' < « < ^/H) 

+ (73) 


The  symbol  «('  labels  the  species  different  from  af . 
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5.  Collision  Term  in  the  Time  Domain 
Klimontovich  and  Puchkin  (197li),  using  a similar  approach 
and  the  same  specialization  to  Eq.  (52),  derive  for  Eq.  (29)  an 
expression  which  they  give  as 

f(»  k-  k'lrjj  Fc 


F = ac'e'n.l-  [A  Ml Zf(Ue) 

-e4-  1 ‘3P  J li*  fr  n 


(7L) 


l€(^k)r  3P 

In  this,  a delta  function  has  been  used  for  the  ion  distribution 
function.  The  variable  is  the  deviation  of  particle  momentum 
from  particle  drift  momentum,  and  the  remaining  time  dependence 
has  been  averaged  over  a period  of  the  driving  field. 
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CHAPTER  III 

HEATING  OF  ELECTRONS  AND  IONS  BY  A HIGH-FREQUENCY  FIELD 

Fluctuat  nec  mergltur. 

Albert  Messiah^ 

In  this  chapter,  the  physical  system  chosen  for  illustration 
is  the  simplified  one  of  a homogeneous  two-component  plasma 
driven  by  an  alternating  electric  field  with  no  magnetic  field. 

It  is  shown  that  the  energy  absorbed  from  the  field  by  one  species 
can  be  expressed  as  the  sum  of  two  terms,  one  representing 
resistive  heating  and  the  other  representing  an  exchange  with  the 
other  species.  The  terms  are  evaluated. 

1.  Introduction 

A hydrodynamic  equation 

The  hydrodynamic  equations  are  the  velocity  moments  of  the  kinetic 
equation  (II. 1).  The  following  exact  expression  is  easily  derived, 
as  shown  by  Chapman  and  Cowling  (1939,  page  49,  Eq.  3. 13, 2): 

+u<h:n7  + nfv-ic  + vn<fiv 

This  equation  has  been  simplified  by  taking  * to  be  a function  of 


Messiah  1962,  front  plate  for  Volume  II. 


L_ 


I 


only.  The  overbar  labels  the  average  with  respect  to  the 
distribution  function. 

.f  j t — P A.  \ 


ij'  Wv)j  - 


I(jf,V»t)  is  the  quantity  which  customarily  appears  on  the  right-hand 
side  of  the  kinetic  equation,  the  "collision  term,"  as  in  (II. 1)  and 
(II. 2)  of  tne  preceding  chapter. 

Deviation  of  velocity  from  drift  velocity: 

V * y - (3) 


Drift  velocity: 


UfSjt)  2 n'1  JSrir 


This  definition  of  u(x,t)  makes  it  the  actual  drift  speed,  which  can 
differ  from  the  u * Jtt  dt  used  elsewhere  in  this  thesis  because 
of  the  collision  term  and  possibly  an  initial  value.  The  difference 
will  be  important  in  a more  general  computation  of  transport  rates; 
in  this  thesis,  the  difference  is  ignored. 

If  there  are  no  space  gradients,  and  if  <f>  is  taken  as  the 
particle  energy  £mV  , then  (1)  reduces  to 

r r i ^ — h o / c \ 


The  separation  of  terms 

In  the  absence  of  the  driving  field,  if  the  system  were  to  begin 
with  nonthermal  distribution  functions,  the  electrons  would  become 
thermal  in  a time  described  by  an  effective  electron-electron 
relaxation  time,  then  the  protons  would  relax  in  a time  greater  by 


Che  root  mass  ratio,  and  finally  Che  temperature  difference  between 
the  two  species  would  relax  in  a time  greater  by  another  root  mass 
ratio  (Spitzer  1962,  p.  136;  Delcioix  1960,  p.  113).  The  interspecies 
equilibration  rate  is  given  by  (1.1)  with  the  relaxation  time 


(111.49): 


Jff,  _ % - & 

c/i  " 71, 


(I.D 


If  a weak  driving  field  is  imposed,  the  species  move  in  opposition, 
and  collisions  convert  some  of  the  energy  of  oscillation  into  thermal 
energy.  Because  each  kind  of  particle  receives  the  same  impulse  in  a 
two-particle  encounter,  the  increment  of  randomized  energy  between  the 
two  species  differs  by  the  mass  ratio,  so  that  the  ion  population  heats 
more  slowly  than  the  electron  population  by  the  mass  ratio.  In  a weak 
field  with  Bt  = the  thermal  energy  of  each  species  grows  at  the  rate 

i***  = <6> 

where  fit  % , the  dissipative  part  of  the  conductivity,  is  given  by 


=-,K,  <’> 

where  V is  the  electron-ion  collision  frequency  (1.8),  and  K is  e or  i, 

O 

The  concepts  of  interspecies  energy  exchange  and  resistive 
heating  can  still  be  given  meaning  in  a system  with  a temperature 
difference  and  a strong  driving  field.  The  rate  at  which  a species 
gains  energy  can  be  derived  from  the  hydrodynamic  equation  (1). 
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If  the  Kilmontovich  collision  term  (II. 6)  is  used  in  (2)  for  the 
left-hand  side  of  (5)»  then  (S)  is 


H 

(8) 

* £•(!£■&)!?> 

(9) 

K 

H 

- <*£■!/„>  -4,-W£> 

(10) 

by  appropriate  definition  of  f j and  Sf  , and  with  species  subscripts 
included  for  the  last  line. 

We  will  find  that  the  second  term  on  the  right-hand  side  is 
quadratic  in  field  strength  for  weak  fields,  so  that  it  can  plausibly 
be  defined  as  the  resistive  heating.  In  the  absence  of  a field,  the 
first  term  is  given  by  (1.1),  so  it  can  be  defined  as  a heat  exchange 
term  These  definitions  aren't  unique;  e.g.,  any  quantity 
(0^  - 0#</)E2f(t)  could  be  added  to  one  term  and  subtracted  from  the 
other  without  destroying  the  properties  which  motivate  those  labels. 

Average  over  a period  of  the  driving  field 

Each  of  the  bracketed  quantities  in  (10)  has  an  expansion  of  the 
form  (11.61).  Taking  the  time  averages  gives 

<u) 

and  the  heating  rate 

T rc  E*  * -H7<  f£fO  « “**<>,  *'“**<>3 

« o» 

This  defines  the  dissipative  part  of  the  conductivity  r -fLC . 
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Matrix  notation  for  the  frequency  dependence 

The  expressions  (11.72)  and  (11.73)  for  the  resistive  and 
equilibration  terms  are  given  a cleaner  appearance  with  more  compact 
notation.  With  arguments  k and  (J  understood,  we  can  write 


Kmn  = 

(13) 

a • L + K*  *■  ■ 

(14) 

N = "ml 

(15) 

Jin.  :nr 

(16) 

i -r  xwfjr/w,)} 

(17) 

l . j‘  = c‘l  l-rjk-J,)} 

(18) 

The  underlining  will  sometimes  be  omitted  from  the  matrix  and  vector 
symbols.  We  define  a shift  operator  by 


S « S' 

n 

so  that  it  has  the  property 


(19) 


(sz) 

ntn 


Clearly  true  is 


(20) 


y$ns  - s <2i> 

- - - n© 

in  which  J*  is  the  Hermitian  conjugate  of  J. 

Remembering  that  the  linear  susceptibility  depends  only  on  the 
magnitude  of  k,  we  get  from  (11.71) 


z 

X 


Jl._  J. 


yti-rn  * 


(22) 


32 


33 


A consequence  of  (23)  and  (II  .71)  is 

r S'  (no  sun).  (28) 

In*  fnn 


2,  Resistive  Heating 


Earlier  work 

The  rate  at  which  a plasma  absorbs  energy  from  a weak  high- 
frequency  electric  field  was  calculated  by  Dawson  and  Oberman 
(1962}  Oberman,  Ron,  and  Dawson  1962;  Dawson  and  Oberman  1963). 
Rand  (196U)  derived  a single-particle  absorptivity  but  didn't 
average  over  the  electron  distribution  function.  Silin  (1965b) 
calculated  the  nonlinear  absorptivity  and  reached  an  expression 
equivalent  to  (111.41).  Raw  and  Salat  (1968)  extended  the  linear 
result  to  stronger  fields  by  taking  more  terms  in  the  expansion  of 
Bessel  functions  with  argutents  proportional  to  field  strength. 
Kidder  (1971)  reaches  a simple  expression  for  the  absorption  of  a 
circularly  polarized  wave.  A quantum-mechanical  derivation  by 
Seely  and  Harris  (1973)  used  an  erroneous  approximation  and  quasi- 
physical arguments  for  the  evaluation  of  an  integral,  and  reached 
a conductivity  which  lacked  the  Coulomb  logarithm  (McKinnis  and 
Goldman  1975).  Klimontovich  and  Puchkin  (197U)  give  a result 


similar  to  that  of  Silin 


Simplified  expression 


Equation  (11.72)  in  the  compact  notation  is,  because  only  the 
real  part  contributes, 

fVyi 

4 5 ^ K S Nj , / (30) 

<w  *'  «r  t ( J 

- i J £ [-ur*'-  (^J^r[h’s  *sk(]nj 

- e‘tf‘  (31) 


Expression  for  unrestricted  k 

This  expression  can  be  developed  further,  but  we  will  use  only 
the  large-k  form.  Silin  (1965b)  and  Klimontovich  and  Puchkov  (1974) 
give  expressions  which  are  similar  to  the  one  which  would  follow 
from  (31). 

Large -wave -number  form 

2 2 

The  susceptibilities  are  products  of  kD/k  with  bounded 
quantities,  so  that  both  M and  N become  unity  at  large  wave  numbers. 
The  fractions  in  (31)  multiply  terms  which  are  of  comparable 
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size,  so  the  ion  term  dominates.  The  contributing  term  is 

lsX'L  - k;0 


(32) 

(33) 


/?<  = -*  /If 

*,/*/  ^ ^ 

1 &'[*/»■*%)-?/**(*.}] 


The  term  in  brackets  is  an  odd  function  of  X , say 

£/  - - r/j;  >s. 

j A-i  A A ^ a A-i  A 

by  changing  the  sign  of  ✓{ , and  also 

--  Z-r  x /) 

^ A Ah  A 

from  the  symmetry  of  the  Bessel  function. 


(34) 

The  sum  is 


(35) 


v -f2T»/££  »r/  ) 


(36) 


~ - t A.  a rz 
' X k-it** 

From  (12),  the  absorptivity  is 

r.  F~"  f ** 

(37) 

The  argument  in  Y^  doesn't  contribute  much  to  the  integral,  so 


C 

b 


■ &(hyj  f <3s> 
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Numerical  value* 

Fig.  3 and  Table  1 show  values  of  Che  dissipadve  pare  of  Che 

ne1 

conductivity  in  units  of  >£ . It  should  be  noticed  Chat  is 0 

itself  depends  on  E through  the  Coulomb  logarithm. 

In  the  measurement  by  Brownell  et  al.  (1974  and  1975),  a value  is 

seen  of  6*(.39)  * 0.42,  differing  from  the  Salat  and  Kaw  value  0.78. 

This  computation  gives  5*(.39)  = 0.87,  which  differs  even  more  from 

o 

Che  experimental  value. 


Fig.  3.  Real  (dissipative)  part  of  the  AC 

conductivity  in  units  of  its  zero-field 


3.  Interspecies  Thermal  Energy  Exchange 
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A 
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Large-wave-number  fora 

The  susceptibilities  Yt(to)  are  (k^k2)  times  quantities 
which  are  bounded  aboTe,  so  that  M and  M are  unit  matrices  at  large 
k. 

<W,.  «, 

Reduction  of  terms: 

Cji KeL  °(l2Ke)  'JZ.K'  * 5 ‘'f'C  <56> 

J ® 

i J?vKnJ7  --  &&£(**%)£*•"*/**!) 


Cj 


- T J*  Jfo)  (57) 

* ° 

*!• 

The  rate  of  heat  transfer  from  ions  to  electrons  is 

<rrt-‘£\,  - 

r-  r!r  v/  , 6L  T»yefci-»K)  (58) 

- > 4/X  i/n/.M  -S 

n i Cv+rUA/* 

= r.  >>; 


and 


Because 


Ke(b/tH(J0) 

Ct>*nu/t 


(59) 

is  an  even  function  of  its  argument, 


is  almost  a delta  function  as  a multiplier  of  it,  we 
Co 

can  put 


c**hu/4 


Co  ~ nut ^ 


*>  (n4/„J 


n<*/_ 


£u;  d — _£<W| 


1 , , 

* [ - "V.tu/c kv-') 


(60) 
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where  the  top  auantity  in  brackets  corresDonris  to  the  too  sign. 

<*;•**  V.  <“> 

Applying  the  k limit  (1.34)  and  then  the  k sum  (1.7)  gives 

<*£•*£>  = f S-Hi <•» 

I 

- Z -55e£_  ±L  V'SMt  i!4j. 

"V  *^cVt  *'  »e'2  e/itfij/l'  '•  "J 

i « 

= l «*> 


Relaxation  rate  in  vanishing  fieldi 


When  E ■ 0,  • I and  r(ft)'°j  so  that 

<‘C  «f>  = ^ 3ny. 


which  agrees  with  (U5). 


Asymptotic  form  in  a strong  field: 


The  asymototic  form  follows  from  (1. 1*0): 


P *v  t(z'tr)  $ 

e 


r ~ *i 

i 

'i  ^ p 


3 


Numerical  values 

For  tabulation,  (64)  can  be  written  as 

<s±-sl>  • s»£y.  {*-%  -Q  Tee)]  <«> 

where  J^(E)  and  J*(E)  ®re  defined  by  equivalence  of  the  expression  to 
(64).  The  equilibration  rate  for  several  values  of  the  electron-ion 
temperature  ratio  is  shown  in  Fig.  4 and  Table  3;  as  for  all  these 
functions,  it  must  be  remembered  that  there  is  an  additional  dependence 
on  field  strength  through  the  Coulomb  logarithm. 


0.0  .5  1.0  1.5  2.0  2.5  3.0  3.5  4.0  4.5  5.0 


Fig.  4.  Electron  heat  loss  rate  in  units  of  3n«n»£  t(&, >» 
for  several  values  of  • 
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Table  2.  Numerical  values  corresponding  to  Fig.  4:  Rate  of  loss  of  electron  heat  to  ions 


CHAPTER  IV 


THERMAL  CONDUCTIVITY 


Ego  in  mi lie  £ormas  transmutavi,  sed  ope ram  meam 
improbum  Problema  perpetuo  luait.  ^ 

James  Bernoulli1 


1.  Introduction 

The  usual  approach  to  calculating  quantities  dependent  on  small 
deviations  of  the  distribution  function  from  equilibria  is  that  of 
Chapman  and  Enskog  (Chapman  and  Cowling  1939;  Grad  1960).  It  Involves 
defining  a smallness  paraiseter 


^ <» 
j expansion  of  the  distribution  functions  in  powers  of yUy 

I 

j 

| and  writing  the  kinetic  equation  in  the  form 

I <2) 

I 

i The  method  is  presented  by  Chapman  and  Cowling  as  an  expansion  to  all 
orders  in  JJL,  The term  requires  Zl^Cf^  .f^)  to  vanish  and  thus 
determines  f<#  to  be  locally  Maxwellian.  The  zero-order  term  is  the 


^Quoted  by  Watson  (1944),  p.  1. 


one  which  describes  Irreversible  processes  In  a near-equilibrium 
plasma  (deGroot  1951),  and  In  practice  is  usually  the  highest  one 


46 


calculated.  The  collision  term  is  bilinear  in  fK  and  tp  , so 

■CJ.'C'V  - r r*  i/'y} *TM, U<3) 

which  is  a linear  integral  equation  in  If  and  subject  to  attack  by 
classical  methods.  The  if  's  are  expressed  as  infinite  sums  of 
orthogonal  polynomials,  and  an  appropriate  integration  gives  a matrix 
equation  which  can  be  inverted  to  a solution  for  the  coefficients  of 
the  polynomials. 

The  thermal  conductivity  to  be  derived  in  this  chapter  can  be 
defined  as  the  K in 

1 « - HCf)-  79.  (4) 

Such  an  expression  can  be  expected  a priori  to  apply  if  jju  is  small, 
as  we  assume  here.  j £ is  the  heat  flux 


fc  = J Jr  ). 


It  is  half  the  quantity  called  the  "thermal  flux"  in  the  book  by 
Kllmontovlch  (1967).  A dimensionless  form  of  S is  convenient: 


- — 

If  the  "Krook  model"  value  (Boyd  and  Sanderson  1969)  of  the  conductivity 
is  used  as  a unit, 


> 


(7) 


then 
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(8) 


end 


<s) 

If  the  only  anisotropies  of  the  plasma  are  the  driving  field  E 
and  the  temperature  gradient,  K will  have  the  form 


<v 


!dt »v 


(10) 


where  and  play  the  roles  suggested  by  their  labels  because  it 
follows  directly  that 

\ = - *'K‘7e  = -Ku(E-*e)  * 

and 

4*  - -k  • ps  . 

Because  the  Interaction  cross-section  for  a two-particle  collision 
decreases  with  the  speed,  and  adding  the  diffusion  velocity  to  the 
driving  velocity  is  more  effective  when  the  two  are  parallel,  we 
expect 


K.  > X. 


a 


In  this  chapter,  a version  of  the  Chapman-Enskog  method  is  used. 
The  thermal  conductivity  is  computed  for  a field  strength  which  is 
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constrained  only  by  nonrelativlstlc  dynamics  and  motionless  Ions.  We 
show  that  the  thermal  conductivity  implied  by  the  collision  term  of 
Klimontovlch  has  the  behavior  as  a function  of  field  strength  which  is 
sketched  in  Figure  S. 

2.  Earlier  Work 

Landshoff  (1949,  1951)  calculated  the  thermal  conductivity  of  an 
equilibrium  plasma,  and  the  other  coefficients  in  Eq.(1.2),  with  a 

fifth-order  Chapman-Enskog  expansion.  Spltzer  et  al.  (Spitzer  and 

I 

j Harm  1953;  Cohen,  Spltzer,  and  Routly  1950)  solved  numerically  a 
second-order  differential  equation  for  a Fokker-Planck  collision  term 
for  an  exact  result.  After  electronic  computers  became  more  useful, 
other  workers  (Hochstlm  1969  and  references  therein)  did  large-order 
matrix  inversions  with  inclusion  of  dependence  on  magnetic  field. 

The  accepted  expression  for  thermal  conductivity,  including 
allowance  for  the  thermoelectric  effect,  is  given  by  Spitzer  (1962): 


Heat  flux  limited  only  by  electron-electron  collisions 


Heat  flux  limited  by  collisions  with 
both  electrons  and  ions. 


Fig.  5.  Sketch  describing  expected  behavior  of  thermal  conductivities 
perpendicular  and  parallel  to  the  driving  field.  The  unit  for  the 
conductivities  Ka  and  is|~S.  The  asymptotes  are  the  values  from 
a 3 x 3 matrix  inversion  known  from  Landshoff  (1949). 
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The  resulc  by  Kelleher  and  Everett  (1967)  that  purports  to  be  the  same 


seems  to  be  in  fact 


ii  « UZI. 


The  value  is  the  one  which  would  result  from  use  of  the  Grad 

t/F 

approximation  for  the  electron-electron  collisions--or,  equivalently, 
a two-dimensional  matrix  inversion  in  the  Chaproan-Enskog  method. 

A recent  measurement  (White,  Kilkenny,  and  Dangor  1974)  involved 
a laser  beam  with  field  strength  which  edged  into  the  strong- fie Id 
regime.  The  result  was 


HBK. 


- O.H  t 0.1. 


This  result  is  in  the  wrong  direction  to  be  explained  by  the 
reduction  of  collision  frequency  by  strong  fields. 

3.  The  Boltzmann  Equation  and  the  Ansatz  for  fe 
Consider  the  Boltzmann  equation 

(ff  +Z-V  - 2-  iyJPcxyrt)  = 

We  can  express  f(x,v,t)  in  the  form 

V re  J 

In  this,  u(t)  is  the  drift  velocity  corresponding  to  the  particle 
acceleration  £(t).  A nondimens ional  velocity  is  used: 

fa)  = e'*** 
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'i  f(v)  and  f0(S)  are  normalized  wich  respect  to  integration  over  their 
| arguments. 

The  quantities  L„(l^)  are  a set  of  orthogonal  polynomials,  chosen 

|j 

i to  be  Sonine  polynomials  (Laguerre  polynomials  of  order  h ).  They  have 


the  property 


| The  first  three  are 


L.W  = I,  L,(t)  -f-t  * F -f*  *t*’  <20) 

The  actual  drift  speed  can  in  principle  be  different  from  u(t). 

The  difference  is  given  by 

i --  P,  w.w.LJMfts) 

* * »/r«  i J * 

• if 

- Z P;  JjtClt)''  ^ iJtJOrrK-1 

= fp( 

- Ci-  <22> 

In  this  chapter,  we  are  computing  the  electron  thermal  conductivity 
with  the  constraint  that  there  is  no  net  current;  which  is  to  say  that 


we  set 


f.  = *. 


By  doing  this,  we  will  reach  a "thermoelectrically  compensated" 
thermal  conductivity  without  an  intermediate  "uncompensated"  value. 

If  there  is  no  net  drift  speed,  then  the  thermal  flux  is  the  same 
whether  it  is  computed  in  the  rest  system  and  then  time -averaged,  or 
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Just  computed  in  the  drift  system.  It  is 


*i  M - s fr,  P"*.  "- 

'Flf, 


r - £ r. . 

2 >* 


(24) 


Thus  the  problem  of  this  chapter  is  to  compute  P^. 

If  f(5f,tf,t)  > f(#,tf-u(t)),  then  the  left-hand  side  of  the  Boltzmann 
equation  (15)  reduces  to  the  gradient  term.  To  first  order  in  ju. , 
only  f contributes.  By  noticing  that 


7 w'  = -w'V  Jn& 


(25) 


we  can  write 


vAa*f?(ziri)  = + VAn,h9  +t*xVAn& 


(26) 


Time  averages  are  taken  by  replacing  v by  v£w  4-  u(t)  and  averaging 
over  the  remaining  time  dependence.  Thus  v-yf  — > ajfcr -rf.  to  lowest 

order  in  and  the  LHS  of  (15)  is 


i 


^ ^ • [ l,  vAk^  -4  rAnfi]  &*')■ 


(27) 


Klimentovich  and  Puchkin  (1974)  applied  such  a time  average  to 
(11.29)  and  give  as  the  result  (11.74),  recopied  here  in  the  form 

l 

i 

v (k  y 


(28) 

(29) 


L 
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4.  Solution  for  the  Heat  Flux 

The  coefficients  of  the  Lj,  on  the  left  of  (27)  can  be  isolated 
by  operating  on  (15)  with  J d^ww^L^ (^w^) . The  result  is 

t j^iA r;  oo> 

The  next  tew  steps  show  how  to  find  the  wanted  coefficient  in 
terms  of  79. 

In  (30),  we  use  (16)  with  7 € { i.l]  only.  The  integration  of 
the  delta  function  in  (30)  is  simple  in  a coordinate  system  chosen  so 
that  k m (0,0, k) . Vector  components  in  this,  the  "internal"  coordinate 
system,  are  related  to  those  in  the  "external"  coordinate  system, 
which  is  fixed  in  space,  by 


zf*  rt..*/ 

1 *j  j 


(31) 


where  ^ has  the  properties 

2 •&«), 


(32) 


and 


A/  * *,-•*/  ^ 


(33) 


The  external  coordinate  system  is  chosen  with  the  third  axis  parallel  to 
the  driving  field.  To  allow  expressing  all  vector  components  in  the 
integrand  of  (30)  with  respect  to  the  internal  coordinate  system,  we  put 
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With  these  change*,  the  first  integral  in  (30)  is 


Definition: 


[$f-  $*- 

»jt  ■ L#]‘  * <37) 

where  H*  is  the  value  contributed  by  Ieft.  This  notation  shortens  (30) 
to 

$,  * £ A OK 


which  is  to  say 


£WF*)  - U*^P  4-  UikP  (39 

t^TTi  * >>  >k  + * ik 

n «.  i/*P  + H*kP  (40 

0 hu  ik  + ik 

The  tensors  can  be  represented  in  the  form  (10): 

= v ® £r  -y.  $7  ♦ r,, , (41 

• jb'V1'  4k  4S*kr  /*v  4*  kS  > 

where  two  ideograms  are  introduced.  In  terms  of  them,  Pj^  is  given  by 

<41‘ 

o =9aP,  <42‘ 

n&i*®upu+&«i,u  (43< 

o -Wu'Vu 
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5.  Evaluation  of  { J ; Explicit  Result  for  the  Conductivity 
| The  electron  contribution: 

Values  for  in  (37)  are  given  by  Landshoff  (1949): 


{ ^ } * K**  [>* 


s*)*  «4* 

Ugorithm  for  evaluation  of  the  Integral: 


Pros  (36), 


[ = / 'a 


The  quantities  (47)  have  the  form 


3 * ftWvfrKWj,*''- 

(j>  (£)  = 

Zr*  Jv*S  L >*S  *Z-VJ 

P.  - ne(z^c:tw% 

V - —n——  L**  T J* -L. 

" m'v*'  " n » 


3 r '£?h« 

Change  of  variables:  i sr^)  jn*  - £ (53^ 

Abbreviation:  V * yn1 

3 ’ £"'§<tr,(!’0'c&k  %A-lJte  m 

Replacing  the  /£  product  by  its  azimuthal  average  allows  use  of 
the  recipe  (1.7)  for  the  part  preceding  the  ft  product^  and  it  becomes 

The  pert  of  (51)  efter  the  ft'*  le.  wlch  f ~enlng  j^jr[\ 

J J'.ll 


I 
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JTj ’* (p  (tv)  = 7L^C'Yrire'T<f* rr»y),  (55) 

n n im  n n *6  2** 

The  effect  of  on  the  polynomial  ^ is  to  convert  T*  to 

" 2 my 

fjfi)  , and  the  effect  ot  HS^C  is  to  convert  Vf  to  . Thus 


7=  Q*»(r$*!r))' 
where  fy&n  is  obtained  from  (f>  by  this  procedure : 

(1)  Take  the  azimuthal  average. 

(2)  Make  the  replacement  wxmzt. 

(3)  Make  the  replacement  t * *r*v» 

(4)  Make  the  replacement  The  result  is  Q(y)  . 

(5)  Make  the  replacement  Yr— > £ to  get  Q(T)» 

The  result  of  applying  the  algorithm  will  have  the  form 


(56) 


(57) 

with  algorithm  applied. 


Because  the  effect  of  an  average  over  the  k-azimuth  is 


K 


we  have  from  (47)  and  (56)  with  (57) 


(58) 


That  is  to  say 


0* 

/*>V 


s 3 


7» 


HV' 


(60) 


1 


(61) 


I 
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Evaluation  of  H In  terms  of  T symbols: 

The  evaluation  of  (56)  will  be  sketched  here.  The  form  for  (57) 

la 

(l^ju  -C'vs; 


B,  ®,  • I*L, 


'A 

cM  = v 


Using  the  expressions  for  and  gives  the  intermediate  result 


I! 

ii 

ii 


0*L,)<i  *£i)J 


(63) 


, , jv  i,i,|  r i 

If**]  ° [l,l,  ifj  ~l  1/1’H  u,)J 

r «i,  111,11, 


(64) 


(64) 


Substitution  of  (20),  the  explicit  forms  for  Ly,  leads  to  the  explicit 
t-dependencei 


{ ‘ 


11  21  12  22 

7/2 

-1 

49/4 

-7 

1 

63/8 

-9/2 

1/2 

441/16 
-189/8 
25/4 
, -1/2 _ 

(65) 


The  first  column  of  (65)  is  to  be  read  = 2v(t-v)(?  - t). 


i i 

I ! 
• 

; 

M 

i 
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I 


2v(t-3v) 


2v(t-v) 


(66) 


Replaclag  t by  T*v  and  then  Tn  by  ni  gives  polynomials  in  v: 

11  21  12  22 

(Z^y)  * 

The  £irst  column  in  (67)  is  to  be  read  « 3 - 2f^,  and  its 

2 

precursor  was  the  polynomial  3v  - 2v  . 


(68) 


Special  values  of  v can  be  used  to  check  the  algebra:  Using  v«0 

in  (65)  and  (66)  gives  the  top  lines  of  (67)  and  (68).  The  top  line 
of  (68)  is  the  same  set  of  numbers  as  Landshoff's  for  zero  field; 
notice  that  a constant  value  of  F will  pass  unchanged  through 
(60)  and  (61). 
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Setting  v **  1 in  (65)  and  (66)  leads  to 


t i 1 *l/4 

fa*'  ,n  i5/8 

b 


* 

% 


9/4  33/16 


77/16  631/64 


correctly  predicting  the  column  sums  in  (67)  and  (68).  Similarly, 
putting  v = -1  in  (65)  and  (66)  leads  to 


* 


33/2  275/8 


fa}  = , 


203/4  1865/16  (?Q) 


2285/16  21563/64 


agreeing  with  the  corresponding  reductions  of  (67)  and  (66)* 

Result  for  the  conductivity: 

The  expressions  for  and  follow  from  comparing  (9)  and  (24) 
with  (44)  and  (45): 

= -1  ^ - (71) 

2 00 

« n t/ 

and  the  corresponding  expression  for  Kjj.  The  quantities  ^ are 
given  by  (59),  (60),  and  (61)  with  (46),  (67),  and  (68). 

Thermal  conductivity  in  zero  field: 

At  E * 0,  the  only  contribution  is  from  fYo)  s K and  (71) 
becomes 

. . s r..s  \ / 67  + !VS  1 Z/r  433  + )8oJT  H “ ' 

Ki  * \ i 

z K#y-/[4.CC¥  -Z.C87]  * A*77  (72) 

= UC  X, 
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Use  of  the  thirteen-moment  representation  of  Grad  for  the 
distribution  function  is  equivalent  to  retaining  only  the  first 
term  in  the  denominator  of  (72),  and  is  thus  unsatisfactory 
for  computing  the  thermal  conductivity. 
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CHAPTER  V 


CONCLUSION 

Ich  glaube,  die  Wellen  verschlingen 
Am  Znde  Schiffer  und  Kahn; 

Und  das  hat  mit  ihrem  Singen 
Die  Lorelei  getan.  , 

Heinrich  Heine 

Three  transport  coefficients  have  been  calculated  as  functions  of 
field  strength  for  a plasma  driven  by  an  electric  field  with  a 
frequency  much  greater  than  the  electron-ion  collision  frequency,  and 
with  a strength  which  does  not  drive  the  electrons  relativistic  but  is 
not  constrained  to  drive  them  at  less  than  the  thermal  speed.  The 
technique  used,  which  involves  application  of  the  Klimontovich  formalism 
to  evaluate  the  collision  term  and  then,  at  an  appropriate  place, 
expanding  an  integrand  in  Inverse  powers  of  the  wave  number,  is  applicable 
to  the  calculation  of  other  transport  properties  of  the  driven  plasma. 

There  are  several  ways  in  which  the  calculation  could  be  generalized. 
Inclusion  of  an  ambient  magnetic  field  is  conceptually  simple  and  might 
be  analytically  tractable.  Other  tra&sport  coefficients,  the  (DC) 
electrical  conductivity  in  particular,  could  be  calculated.  Allowing 
a more  general  polarization  for  the  electric  field  is  probably  simple 
and  direct;  in  particular,  the  response  to  the  electric  field  is  simpler 

■ 1 - 11 


^From  "Lorelei." 


for  a circular  polarization  than  for  linear  polarization  because  the 
driven  speed  of  an  electron  is  constant.  Including  relativistic  effects 
by  perturbation  or  otherwise,  would  extend  the  permitted  field  strength. 
The  use  of  a quantum-mechanical  kinetic  equation  is  an  intriguing 
concept;  it  is  known  to  be  feasible  for  calculating  inverse  Brems- 
strahlung.  A quantum-mechanical  formulation  might  allow  the  application 
to  more  general  statistics.  The  calculation  of  the  second  terms,  of 
order  unity,  in  the  several  expansions  which  begin  with  In  A i and,  in 
particular,  including  powers  of  k past  the  first  in  the  expansions 
beginning  k-1(k“^  +. . .);  would  give  a quantitative  estimate  of  the 
limits  of  validity  of  the  results. 

The  polishing  and  re-organizing  of  the  problem  for  simplest 
presentation  and  greatest  usefulness  of  the  results  is  much  less 
complete  than  the  author  would  like.  All  of  the  topics  mentioned 
above  as  possible  extensions  might  have  been  followed  to  at  least  the 
first  calculational  difficulty,  with  an  expectation  of  including  most 
of  them.  To  give  a general  expression  for  integrals  of  the  collision 
term  and  express  the  particular  results  as  special  cases  would  have 
been  more  pleasing.  In  the  results  for  the  thermal  conductivity, 
there  appear  several  linear  combinations  of  the  quantities  Pr(E); 
most  of  these  combinations  have  the  property  that  the  asymptotic 
behavior  is  E"^  even  though  the  individual  P's  go  to  ONE"* + . . . 
because  the  leading  terms  cancel.  This  behavior  suggests  the  existence 
of  an  organizing  principle  not  recognized  in  the  calculation.  The 
expression  for  the  inverse  Bremsstrahlung  rate  can  be  reduced  to  a 
reasonably  simple  form  before  use  of  the  large-k  limit;  it  was  omittei 
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because  the  result  isn't  new.  A better  overview  of  the  Klimontovich 
formulation  and  a comparison  with,  e.g.,  the  BBGKY  hierarchy,  could 
have  been  included.  Is  there  a paradox  lurking  in  the  frequency- ratio 
factor  in  the  Coulomb  logarithm,  and  how  is  the  transition  made  to 
the  vanishing- fie Id  value?  Does  the  frequency  ratio  explain  the 
discordant  result  of  White,  Kilkenny,  and  Dangor? 


CHAPTER  VI 


LIST  OF  SYMBOLS 

When  a page  number  is  given,  it  locates  the  definition  in  the 
test,  or  shows  where  the  symbol  is  used.  Most  of  the  various  integer 
indices  and  dummy  variables  are  not  mentioned. 


Non- alphabetical 

{ y Ensemble  average,  p.  14. 

k,  k These  two  labels  for  a vector  are  used  interchangeably. 

A (1)  Average  with  respect  to  distribution  function. 

(2)  Average  of  A over  a period  of  the  driving  field. 

( } A label  for  the  part  of  — ; — attributable  to  noncollective 

' 3t  Lll  Jt 

interparticle  interaction,  p.  13. 

M Hermitian  conjugate  of  the  matrix  M. 

Components  of  the  matrix  (jf/  UQ),  p.  54. 


i 

i 
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Roman 

Hit 

& (IX  t)  "Microscopic"  acceleration.  The  exact  acceleration 

of  a particle  with  coordinates  (x,  v,  t)  in  a particular 
member  of  the  ensemble.  Used  p.  14. 

r 

(X>  Ensemble  average  of  £ , p.  14. 

A particle  acceleration  corresponding  to  imposed  force. 

Used  p.  14. 

bQ  Classical  distance  of  closest  approach,  p.  7. 

c Speed  of  light. 

r«f, 

d^  Excursion  distance,  d„  mff^l » P*  22. 

E (1)  Nondimensionalized  field  strength,  p.  11. 

(2)  An  unspecified  field-strength  variable. 

E0  Amplitude  of  the  sinusoidal  electric  field,  p.  21. 

E^^AE),  y<E)  Functions  which  appear  in  the  integrands  of 

expressions  leading  to  thermal  conductivity,  pp.  57  and  59. 
e Proton  charge . 

"Microscopic"  distribution  function.  The  exact  distribution 
function  for  a particular  set  of  point  particles,  as 
distinguished  from  the  ensemble  average.  P.  13. 
f (xvt)  Distribution  function. 

£(t)  Arbitrary  or  unspecified  function. 

G(x,x' »v,v' ,t,t')  Part  of  a Green's  function,  p.  15. 

Gj^VjV' ,t,t')  Part  of  a Green's  function,  p.  16. 
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! 

I 


I (vt)  or  I(vt)  Collision  term  in  the  kinetic  equation  for  species *( 
or  for  unspecified  species,  pp.  13  and  17. 

I (1)  The  energy  flux  of  an  electromagnetic  wave. 

(2)  Collision  term. 

I >J  Momentum  integral  of  ^(U)  weighted  by  the  collision  term, 

p.  28  • 

In(x)  Bessel  function,  used  p.  10. 

Jn(x)  Bessel  function. 

J Electrical  current,  p.  3. 

J,  Je  Column  vector  with  elements  P*  31. 

J^,  L Column  vector  with  elements  JnQ£*d^),  p.  31. 

K (1)  Thermal  conductivity,  p.  3. 

(2)  Thermoelectrically  compensated  thermal  conductivity, 
p.  49  ff. 

K0  Thermal  conductivity,  reference  value,  p.  46. 

K^,  Components  of  thermal  conductivity  perpendicular  and 

parallel  to  the  field,  p.  47. 

K Susceptibility,  p.  34. 

k Fourier  transform  variable,  p.  19. 


The  "Debye  wave  number"  • k 


i ^ 


. /-fsajti ) 

Otf  ' & ' * 

• tt 


^max*  ®Snin  Empirical  limits  on  k,  p.  5. 

9(  rs  Weight  function  in  the  collision  integral,  p.  52. 


M «« 
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l;\x) 


M 

M(a,  b. 


M 


m 


n,  n„ 


<? 


s, 

A 

S^Cxvc) 


Liouville  operator,  p.  16. 

Sonine  polynomials,  p.  51. 

Dielectric  matrix,  p.  29  and  34. 

-2i^)  Confluent  hype rgeorae trie  function,  p.  10. 

Same  as  m^. 

Mass  of  a particle  of  unspecified  species,  or  of  an  electron. 
Mass  of  a particle  of  species  * . 

Inverse  of  dielectric  matrix,  p.  31. 

Number  density  for  unspecified  species  or  electrons  or 
species  oc  . 

The  "order"  symbol.  y •&  x means  that  lul^  /?"/ 
is  finite. 

Momentum  corresponding  to  the  deviation  velocity  U. 

Pressure  of  one  species,  P * nd. 

Stress  tensor.  P^j  » Jj'v’  P. 

Traceless  stress  tensor. 

Coefficients  in  an  expansion  in  Sonine  polynomials,  p.  51. 
Rotation  operator,  p.  53. 

Thermal  flux,  p.  3. 

Dimensionless  thermal  flux,  p.  46;  components,  p.  47. 

Source  term,  pp.  15  and  17. 

Shift  operator,  p.  34. 

Same  as  S. 

An  initial  time.  It  precedes  greatly  any  other  times 
defined  by  the  stationary  system,  but  is  recent  enough 
that  collisions  have  not  changed  the  distribution 


function 
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u 

Velocity  relative  to  drift,  p.28  . 

Amplitude  of  driven  velocity,  p.  22. 

uQc.t) 

Drift  velocity,  p.  28. 

Veff 

P.  7. 

V 

Velocity. 

ve»  v 

Thermal  speed,  v^  « 

VD»  VE 

V " r 

Electron  driven  speed,  vE  * ■■■  ~ . 

w 

Dimensionless  velocity,  p.  50. 

x°Qc,  v,  t > 

Orbit  coordinate,  p.  15. 

2Si(t) 

Space  coordinate  of  the  i^  particle 

Z 

Charge  number  for  sP®cies  * • z±m 

Z 

Ion  charge  number. 

p-13  • 


1. 
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Greek 

oc  Fine  ^ructure  constant,  used  p.  7. 

*-tP  (1)  Species  subscript,  e or  i. 

(2)  A thermoelectric  coefficient,  p.  3. 
e('  The  species  different  from  Of . 

r/s)  A function  which  appears  in  the  analysis  in  Chapters 
III  and  IV,  pp.  8 and  10. 

A^ tt  * Velocity  acquired  between  times  t’  and  t,  p.  23. 

$ Operator  on  a field  quantity  which  defines  a fluctuation, 


p.  I]*. 


if 

€ 

en 

47 

®(t) 


Defined  on  p.  U|if  each  letter  f is  replaced  by  2.  . 

An  electric  field  related  to  fQr  by  (xvt)=-  — jE(xvt). 

*V  " 

P.  17. 

Jjyj vf  r , p . 30 . 

P.  18. 

P.  18. 

Watson's  notation;  e = 2 - S\  n,  used  only  in  Appendix  A. 

n n,u 

Convenience  variables,  p.  11. 

The  Heaviside  step  function,  zero  for  negative  arguments 
and  unity  for  positive  arguments. 

Temperature  in  energy  units. 

Argument  of  Coulomb  logarithm,  p.  5. 

Bessel  fungtlpn,  p.  10. 

The  Debye  shielding  distance.  See  kQ. 

de  Broglie  wavelength  of  a thermal  electron,  p.  8. 
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M'  (1)  Dimensionless  parameter  proportional  to  the  gradient, 

p.  45. 

(2)  A spherical  coordinate,  cosine  of  the  colatitude. 

(3)  In  Chapter  IV,  an  integer  subscript. 

^ Electron- ion  collision  frequency,  p.  7.  Often  in  the 

literature  called  "ei- 

1^  Effective  collision  frequency  for  heat  exchange,  p.  43. 

S’  (1)  A thermoelectric  coefficient,  p.  3. 

(2)  AC  electrical  conductivity,  p.  29. 

Dissipative  (real)  part  of  f?  , p.  30. 

T ^ Electron-ion  equilibration  time,  p.  3. 

X(S)  A parameter  .in  the  heat-exchange  rate,  p.  43. 


y m 


co  u> 

?!  P+ 


Perturbation  of  distribution  function,  p.  46. 

A phase,  p.  17. 

Arbitrary  function  of  U,  p.  27. 

Linear  susceptibility  of  species  of,  pp.2li  and  11. 
Susceptibility  contributed  by  species  Pi,  p.  18.  There  are 
two  susceptibilities,  distinguished  by  the  number  of  time 
or  frequency  arguments. 

Frequency  matrix,  p.  31. 

Angular  frequency,  a Fourier  transform  variable. 

Angular  frequency  of  the  driving  field. 

Plasma  frequency  for  electrons  or  species  Of  t = 5-*- 
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APPENDICES 


Appendix  As  Alternative  Derivation  for  Tr(i  ) 
This  appendix  shows  another  derivation  of  (1.26)  for  fj.* 


Lemmas 


ZZ  = Tjjjjr*”  * O'*"  ' <*-l> 


Proofs 


From  Watson  (19l*U,  p.  36),  we  have  Neumann's  expansion: 


(i*r=  0 


(A.  2) 


= £?f  £ 


(mi  ' iVr 

= rtf  ZT  ^rO"')ar-0nir-\.  Jjlfa.)) 


Thus  each  sum  is  a polynomial  of  order  r in  x2  „ The 

highest-order  term  is  x2r  lead8  (A«l)» 


The  definition  (1.10)  gives 


r (i) = imf  (r'Sf.  4V//V^A  (a.u) 

7**#  »>*-•  /».-<>  r*  v / 


Using  the  lemma  makes  this 


I MUO  | 

n /i)  ^fr 


which  is  the  same  as  (1.26). 
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Appendix  B:  Pj.  in  terms  of  Bessel  Functions 

There  is  a formula  in  Abramowitz  and  Stegun  (196k,  p.  506) 
which  was  incorrect  in  early  printings: 


In  this. 


(B.l) 


_ fn+fe-i)/ 


r\  \ - 

‘H  = Ct  -Of  ’ 


(B.2) 


We  put  a ■ q ♦ £,  b ■ q + d ♦ 1,  expecting  to  take  the  limit  as 


d goes  to  zero: 


(B-3> 

The  factors  related  to  the  limit  are  underlined.  If  n ■ 0,  they  become 

rcdU  = ci!  ->  i, 

and  if  n f 0, 

rrj)  h 

(Zcf-I)l  (ef-f-i)!  (j+J+n)! 

= h^L  — A,.n/  *Q...  (a f-9+H’O  [! 

c/ 

_ o„i  iL  vr'J-ft*'-) 

M(rt-H'i-u')- e-**f ♦£ 

= r m"/\  fi , — tt-  (b.i*) 

o (!-*)!  (f»)! 
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A simDle  check  of  that  is  available:  Putting  q n:akes  it 


/I"^V 


I (-/A.KV  - £/u^  (B-5) 


OP 


1 • 1. 

Also  in  agreement  is  K(£,l,-2i2)  ■ Ac^  )• 


Combining  (B.U)  with  the  othtr  factors  of  (1.26)  gives. 

r r*  - ^ *,r 

>• 


o-»;! 


Appendix  C:  Recursion  Formulas 

This  appendix  shows  how  recursion  formulas  from  Abramowitz 
and  Stegun  (1961*  ) translate  to  formulas  for  the  V symbols.  The 
arguments  are  related  by  * ■ -2^2. 

Recursion: 

Setting 

a » b * rt  I 

and  multiplying  by  Sivea 

(i-vrjs)  rj*)  * 5 0 


r-r- 
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Differential  recursion: 


M'^Z ) = (<*+•; 

[x'MC^l)]'*  ^zrM(a*ijii-iJz)  +rz'r',M(a\z) 

o-f)! 

Including  a factor  — t ,y.  '“7  c gives 

I".  (*f7. 

*£<;  - = -rr/«  ^a) 

because  of  Z 1 *2i. 

Appendix  D:  Initial-Time  Correlation  Function 

If  the  two-particle  correlation  function  is  ignored,  Fq.  (U2) 
from  Wu  (1966)  is 

\ * p - * y < — * ~~ 

Taking  the  Fourier  transform  with  respect  to  each  of  x and  x 1 gives 

Ks^/rOrP  ^isV)}  - f ^ S(r-r') pxjVc f frr  t ) 

= tyni(r(r-r'J  rcyt) 

in  which  F (x,v,t)  is  assumed  independent  of  x,  and  the  unit  box 
normalization  has  been  used.  This  is  (11.44). 


ildb 
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Appendix  E:  A Property  of  the  Green 'a  Function 


Claim? 


Q (rv'HI  = JjV  GJntn  G'+CvVi't,) 

* * * 


Proof? 

The  Integral  la 


Using  the  delta  functiona  gives  aa  the  su«  of  the  phases 

* Vi  W'^'t-XO 
; h-  { y'ft-ri  * Jjr*.Tt.}  t 

* i*  { v'(t-t')  -(t-r)A^  ttt  +jf/r*rt,} 

» A-  {v'ft-ty  t jOrd  ri '} 

* *4<XV> 

jrilmi^tAng  v»  • from  the  delta  functions  gives  the  correct  argument 
for  the  remaining  delta  function,  completing  the  proof. 
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Appendix  F:  Rationale  for  Che  Choice  of  k^v 

A crude  argument  is  given  in  two  steps: 

Step  One:  The  kinetic  equation  for  the  fluctuations  is 

simplified  by  ignoring  terms  quadratic  in  the  fluctuations.  These 
terms  are  important  when  particles  are  close.  It  is  known  (vide,  e.g., 
Krall  and  Trivelpiece  1973)  that  close  collisions  contribute  only  a 
fraction ) to  the  Rutherford  scattering  cross-section.  Ignoring 
collisions  which  produce  deflections  greater  than  ^ preserves 
the  validity  of  the  kinetic  equation  without  losing  important  effects. 

Step  Two:  Those  parts  of  a configuration  of  particles  which 

involve  particles  closer  than  a distance  b are  described,  when  the 
space  variable  is  replaced  by  a Fourier-transform  variable  k,  by  values 
of  k satisfying 


so  the  region  correctly  described  by  the  linearized  kinetic  equation 


is  the  one 


